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Introduction: 



The present work includes some of the author’s original researches on the integer 
solutions of equations and linear systems: 

1 . The notion of “general integer solution” of a linear equation with two unknowns 
is extended to linear equations with n unknowns and then, to linear systems. 

2. The properties of the general integer solution are determined (both of a linear 
equation and of a linear system). 

3. Seven original integer algorithms (two for linear equations and five for linear 
systems) are presented. The algorithms are carefully demonstrated and an 
example for each of them is presented. These algorithms can be easily introduced 
into computer. 
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INTEGER SOLUTIONS OF LINEAR EQUATIONS 

Definitions and properties of the integer solutions of linear equations. 

Consider the following linear equation: 

n 

(i) Yj a i x i =b ’ 

i = 1 

with all 0 and b in Z . 

Again, let h eN , and yi : Z* — » Z , i = \,n . ( 1 ,n means: all integers from 1 to n ). 

Definition 1. 

x ; = a ; ° , i = \,n , is a particular integer solution of equation (1), if all a, 0 e Z and 

= b . 

i = 1 

Definition 2. 

a, = f j (k l ,...,k h ) , i = 1 , n , is the general integer solution of equation (1) if: 

a ) i i a i f i (k l ,...,k h ) = b; V(i, Z\ 

1 = 1 

b) For any particular integer solution of equation (1 ),a, = a f ,i = 1 ,n, there exist 

(k° ,...,k°) e Z /! such that a ° = f(k° ,...,k°) for all i = l,n {i. e. any particular integer 

solution can be extracted from the general integer solution by parameterization} . 

We will further see that the general integer solution can be expressed by linear 
functions. 

h 

For 1 < i < n we consider the functions f] = ^ c i) k j + d i with all c jf , d i eZ. 

j = i 



Definition 3. 

A = (c l] ) l j is the matrix associated with the general solution of equation (1). 

Definition 4. 

The integers k v ...,k s , 1 < s <h are independent if all the corresponding column 
vectors of matrix A are linearly independent. 

Definition 5. 

An integer solution is s -times undetermined if the maximal number of 
independent parameters is 5 . 

Theorem 1. The general integer solution of equation (1) is (n - l)-times 
undetermined. 
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Proof: 

We suppose that the particular integer solution is of the form: 

r 

(2) x j = Yj u ie P e + i = hn, with all u ie , v, eZ, 

e=l 

P e are parameters of Z , while a < r < n - 1 . 

Let (x[ x”) be a general integer solution of equation (1) (we are not interested in 
the case when the equation does not have an integer solution). The solution: 

x j =a n k j +x°j, j = l,n- 1 

( " _1 ^ 

X n=-\ll a j k j- X nJ 

is undetermined (n -l)-times (it can be easily checked that the order of the associated 
matrix is n — 1). Hence, there are n — 1 undetennined solutions. Let’s consider, in the 
general case, a solution be undetennined ( n - l)-times: 

n — 1 

x, = y.c.f, + dj , i = 1 ,n with all c ijt d, e Z . 

7 = 1 

Consider the case when b = 0 . 

Then 

n 

2>, = °. 

i—i 

It follows: 

n n ^ n - 1 ^ n n - 1 n 

Zcx, = I ZcA + 4 I = = 0 • 

1=1 i=l V 7=1 / i‘=l 7=1 i - 1 

n 

For kj = 0 , j = 1,77 — 1 it follows that = 0 . 

i=l 

n 

For k : =1 and k, =0, j * j n , it follows that V = 0 . 

7o J J J u * */o 

1=1 

Let’s consider the homogenous linear system of n equations with n unknowns: 

n 

ZlC/ = 0, 7 = !.«-! 

1 = 1 
n 

T J x l d, = 0 
/= 1 

which, obviously, has the solution x ( . = a j , i = l,n different from the trivial one. Hence 
the determinant of the system is zero, i.e., the vectors c j = (c yj ,...,c nj ^ , j = l,n-l, 
D = (d l ,... , d n )' are linearly dependent. 
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But the solution being (n - 1) -times undetennined it shows that c., j = 1 ,n - 1 are 
linearly independent. Then (c 1 ,...,c n _ 1 ) determines a free sub-module Z of order n -1 in 
Z n of solutions for the given equation. 

Let’s see what can we obtain from (2). We have: 

n n f r ^ 

0 = Z a i X i =E fl i I U ‘e P e + V i • 

As above, we obtain: 

n r 

Tj a i v i = ° and Yu a < u < = 0 

i = 1 <?=1 

similarly, the vectors U h = (u Ui ,...,u nh ) are linearly independent, /z = l,r , U h , h = 1 ,r are 
y = (v 1 ,...,v„) particular integer solutions of the homogenous linear equation. 

Sub-case (al) 

U,h = l,r are linearly dependent. This gives {£/,,...,£/ }= the free sub-module of 

order r in Z" of solutions of the equation. Hence, there are solutions from jV' t ,...,L i ,} 

which are not from r ] ; this contradicts the fact that (2) is the general integer 

solution. 

Sub-case (a2) 

U h , h = l,r, V are linearly independent. Then [U 1 ,...,U r }+V is a linear variety 
of the dimension <n — 1 = dim[\/,...,L ,} and the conclusion can be similarly drawn. 



Then: 



n 

Consider the case when b * 0 . So, '^ j a i x j = b . 

1 = 1 



n ^ n - 1 ^ n— 1 f n A n 

Zc I H c u k j + d i I = Z Z a i c n k i + Z a A = b > v<,K-X_i ) e • 



i = 1 V;=1 



j=l v i= 1 



As in the previous case, we obtain = b and Z^ = 0, V j = 1, n - 1 . 

/= 1 i=l 

The vectors c ; = [c ip ...,c nj j , j = l,n - 1, are linearly independent because the solution is 
undetermined (n - l) -times. 

Conversely, if c 1 ,...,c n _ 1 , D (where D = (d l ,...,d n J ) were linearly dependent, it 



n — 1 

would mean that D =/ s c with all s scalar; it would also mean that 

i—t j j j ’ 

j = i 




This is impossible. 



a i c ti 



= 0 . 



(3) Then \c l ,...,c n _ l ^ + D is a linear variety. 
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Let us see what we can obtain from (2). We have: 

n n ( r \ r f n \ n 

b = =Y J a ,{ Ys U ie P e + L I = II Tj a i U ie P e + Z fl / V f 

1=1 1=1 V e= i / e= \ V i = i J ,=i 



and, similarly: 2_ i a l v i = b and ^a ; M ie = 0, Ve = l,r, respectively. The vectors 

/= 1 7 = 1 

U e = (« le ,...,« n( ,y , f = 1 , r are linearly independent because the solution is undetermined 
r -times. 

A procedure like that applied in (3) shows that U v ...,U r , V are linearly 
independent, where V = (v l ,...,v n > f . Then [U l ,...,U r }+V = a linear variety = free sub- 
module of order r<n- 1. That is, we can find vectors from \c 1 ,...,c n _^ + D which are 
not from \U 1 ,...,U r } + V , contradicting the “general” characteristic of the integer number 
solution. Hence, the general integer solution is undetermined (n - l)-times. 



Theorem 2. The general integer solution of the homogeneous linear equation 

n 

^ a l x i = 0 (all cl eZ\{0}) can be written under the form: 

1 = 1 

71 - 1 

(4) x , = H c u k i’ i = l ’ n 

7=1 

(with dy = ... = d n =0). 

Definition 6. This is called the standard fonn of the general integer solution of a 
homogeneous linear equation. 

Proof: 

We consider the general integer solution under the form: 

72 - 1 

x i=H c ij p j +d i’ / = 1 ’« 

7=1 

with not all <r/ ; = 0 . We’ll show that it can be written under the fonn (4). The 
homogeneous equation has the trivial solution x t = 0, i = 1, n . There is 

71-1 

) e Z" such that Yep 0 , + d i = 0, Vi = 1 , n . 

7 = 1 

Substituting: P. = kj + Pj, j = l, n — l in the form shown at the beginning of the 
demonstration, we will obtain form (4). We have to mention that the substitution does not 
diminish the degree of generality as P. e Z <=> kj eZ because 7 = 1,71-1. 



Theorem 3. The general integer solution of a non-homogeneous linear equation is 
equal to the general integer solution of its associated homogeneous linear equation plus 
any particular integer solution of the non-homogeneous linear equation. 

Proof: 
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n — 1 

Let’s consider that x i = / j c jj k j , i = \,n, is the general integer solution of the 

j= i 

associated homogeneous linear equation and, again, let x t = v,, i = 1 ,n , be a particular 

n — 1 

integer solution of the non-homogeneous linear equation. Then x i = H c ij k j +v i ’ i = l ’ n ’ 

M 

is the general integer solution of the non-homogeneous linear equation. 

n n f n - 1 n f »-l A n 

Actually, £ a,x, = £ a,. £ c, ; L ; + v,. = Z a , Z c a k j + Z a » v i = b ; 

i=i i=i V J ,=1 V t=i y i'=i 

if x ; =x°, i = l,n, is a particular integer solution of the non-homogeneous linear 
equation, then x t = x t - v,, i = l,n , is a particular integer solution of the homogeneous 
linear equation: hence, there is e Z" 1 such that 

Z C '<A° = ~ v P V i = 1 ,n , 

j= i 

i.e.: 

«- 1 

Z^ 0 + v > = *° ’ V i = ’ 

t=i 

which was to be proven. 



Theorem 4. If x f = ^ c (> kj, i = \,n is the general integer solution of a 

7=1 

homogeneous linear equation ( c tj , . . . , c, y . ) ~ 1 V j = 1, n - 1 . 

The demonstration is done by reduction ad absurdum. If 3j 0 , 1< j 0 <n- 1 such 
that (c,.. ,...,c . ) ~ d . ^ ±1 , then c,.. = c’., d, with I c,’, . ) ~ 1, V / = 1,/z . 

But x i =c' ij , i = l,n, represents a particular integer solution as 



Z a Ai = Z fl TV„ = 1 ' d j 0 • Z«.- c So = 0 

/=! i=l i=l 



(because x f = c tj , i = \,n is a particular integer solution from the general integer 
solution by introducing kj = 1 and k f = 0 , j ^ /„ . But the particular integer solution 



x,. = c'y , i = l,n, cannot be obtained by introducing integer number parameters (as it 

should) from the general integer solution, as from the linear system of n equations and 
n — 1 unknowns, which is compatible. We obtain: 

n 

x ; = / c-.k.+c,- d : k ; =c .. , i = \,n. 

I i—i V 7 'Jo Jo Jo 'Jo ’ ’ 

7=1 
7 - Jo 

Leaving aside the last equation - which is a linear combination of other n — 1 
equations - a Kramerian system is obtained, as follows: 
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C, , C C, , 

11 I/O 1 



, C n-l,l" ,C n-l /„ ••• c n-ln-l 1 _ 

k : = 7 ; T = 



c, c„ d , c, „ , 

11 Vo Jo 



c , , ...c , d ...c , , 

n- 1,1 n-ljo Jo n-ln-1 



d. 



Therefore the assumption is false (end of demonstration). 

Theorem 5. Considering the equation (1) with (a 15 ...,a n ) ~ 1, b = 0 and the 

ft- 1 

general integer solution x j = i = l,n , then 

j = i 

(a 1 ,...,a M ,a i+1 ,...,a fi ) ~ (c,.,,. Vi = l,n . 

Proof: 

The demonstration is done by double divisibility. 

ft — 1 

Let’s consider *'o> arbitrary but fixed, vy = yV, ,fc, . Consider the 

7=1 

equation ^ a i x i = -a ig x t . We have shown that x ; = cy, i = l,n is a particular integer 



solution irrespective of j, a < j < n - 1 . 

The equation ^ a j x l = ~a ig c igj obviously, has the integer solution x t = cy, i * i 0 . 

ijii 0 

Then [a l ,...,a ig _ x ,a ig+l ,...,a n ^ divides - a jg c igj as we have assumed, it follows that 
(uj,...,a n ) ~ 1, and it follows that (a x ,...,a jg _ l ,a ig+x ,...,a\\c igj irrespective of j. Hence 

[a l ,...,a jg _ l ,a ig+l ,...,a\\[c igX ,...,c ign _ l ), Vi' = l,n, and the divisibility in one sense was 
proven. 

Inverse divisibility: 

Let us suppose the contrary and consider that 3/, e\,n for which 
(a l ,...,a h _ x ,a k+x ,...,a n )~d kl *d h2 ~(cy c.^); we have considered d hl and d l<2 
without restricting the generality. d^\d^ 2 according to the first part of the 
demonstration. Hence, 3d eZ such that d i 2 = d ■ d t , , \d \ ^ 1 . 

ft — 1 ft — 1 

X ■ = / C; f ;=d ■ d; . / C: -k , I 

<1 'i7 J MZ-7 'i7 J ’ 



7 = 1 



7 = 1 



ZcL = 0 => ZcL = 



(, J • d: , 7 C, :k; 
h h 1 h J J 



i=l 



7 = 1 
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where ~ 1 . 

The non-homogeneous linear equation c/ (| , has the integer solution 

i ?- /j 

because a, d. , is divisible by (a 1 ,...,a._ 1 ,a, i+1 ,...,a n ) . Let’s consider that x t =x°, i ^ , 

n 

is its particular integer solution. It follows that the equation '^ J ci j x j = 0 the particular 



i = 1 



solution v,. = x°, i x h = d t , which is written as (5). We’ll show that (5) cannot be 
obtained from the general solution by integer number parameters: 



( 6 ) 



n - 1 



Tj c ij k j =x % 



7= 1 



n — 1 



cl ■ cl c„L ; = cl , 
'i 1 i—. y j h 1 



7 = 1 



But the equation (6) does not have an integer solution because d ■ d t , I , thus, 

contradicting, the “general” characteristic of the integer solution. 

As a conclusion we can write: 

n 

Theorem 6. Let’s consider the homogeneous linear equation = 0 , with all 

i = 1 

a t e Z\{0| and (a p ...,a n ) ~ 1 . 

h 

Let = '^ j c jj k j , i = \,n, with all c tj eZ, all k j integer parameters and let’s 

7=1 

consider h e N be a general integer solution of the equation. Then, 

1) the solution is undetermined (n - l) -times; 

2) V j = 1, n - 1 we have (c 1; , ..., c nj ) ~ 1 ; 

3) V i = l,n we have (c a ,...,c in _ x ) ~ (a 1 ,...,a M ,a i+1 ,...,a„) . 

The proof results from theorems 1,4 and 5. 

Note 1. The only equation of the form (1) that is undetermined n -times is the 
trivial equation 0 • x l + ... + 0 • x n = 0 . 

Note 2. The converse of theorem 6 is not true. 

Counterexample: 



(V) 



X j — — fcj + k 2 

x 2 = 5fcj + 3k*, 



[x 3 = lk x - k 2 ; ipfeZ 

is not the general integer solution of the equation 
(8) -\3x x + 3 x 2 - 4x 3 = 0 
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although the solution (7) verifies the points 1), 2) and 3) of theorem 6. (1, 7, 2) is the 
particular integer solution of (8) but cannot be obtained by introducing integer number 
parameters in (7) because from 

— /Cj + k 2 = 1 

5£j + 3k 2 = 7 
7fcj - k 2 = 2 

1 3 

it follows that & = — g Z and k = — <£ Z (unique roots). 
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AN INTEGER NUMBER ALGORITHM TO SOLVE LINEAR 

EQUATIONS 



An algorithm is given that ascertains whether a linear equation has integer number 
solutions or not; if it does, the general integer solution is detennined. 

Input 

A linear equation a l x l + ... + a n x n = b , with a t ,b e Z , x t being integer number 
unknowns, i = 1, n , and not all a i = 0 . 

Output 

Decision on the integer solution of this equation; and if the equation has solutions 
in Z , its general solution is obtained. 

Method 

Step 1. Calculate d = ( y a l ,...,a n ) . 

Step 2. If d lb then “the equation has integer solution”; go on to Step 3. If d\b 
then “the equation does not have integer solution”; stop. 

Step 3. Consider h := 1 . If \d\ ^ 1 , divide the equation by d; consider 

n, := a i Id, i = 1, n, b :=b I d . 

Step 4. Calculate a = min In I and detennine an i such that a, = a . 

a s *0 

Step 5. If a ^ 1 then go to Step 7. 

Step 6. If a = 1 , then: 

(A) Xj = -(a,*! + ... + + a M x M + ... + a n x n - b) ■ n,. 

(B) Substitute the value of x t in the values of the other determined 
unknowns. 

(C) Substitute integer number parameters for all the variables of the 

unknown values in the right term: k v k 2 ,...,k n _ 2 , and k n _ 1 

respectively. 

(D) Write, for your records, the general solution thus detennined; stop. 
Slep7. Write down all a p j ^ i and under the form: 

a t = c ¥lj + 0 





a, 




b 


b = a cj + r where r/ ( = 


i_ 


, 4 = 






_ a i! 


_ a i_ 



Step 8. Write v, = -cpx x - ... - r/,- .v, , - q i+1 x i+1 - ... - q n x„ +q- t h . Substitute the 
value of x t in the values of the other detennined unknowns. 

Step 9. Consider 
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a , , := r. , 6 := r 

< and 

«,- + i : = 'Ll : = h 

: \ h := h + 1 

«» := r « 

and go back to Step 4. 



Lemma 1. The previous algorithm is finite. 

Proof: 

Let’s fl^ + ... + n„x n = b be the initial linear equation, with not all a, = 0 ; check 
for min In I = a, ^1 (if not, it is renumbered). Following the algorithm, once we pass 

a s *0 11 ^ *“ 

from this initial equation to a new equation: a]x t +a 2 x 2 + ... + a n x n =/?', with |nj| < |n,.| for 
i = 2 ,n , \b'\ < \b\ and a l = -n, . 

It follows that min In < min In I . We continue similarly and after a finite number 

a' s *0 1 1 a s *\ 1 S| 

of steps we obtain, at Step 4, n := 1 (the actual n is always smaller than the previous n , 
according to the previous note) and in this case the algorithm terminates. 



Lemma 2. Let the linear equation be: 

(25) n.x, + a 2 x 2 + ... + a x = b , with min In I = n. and the equation 

a s *0 

(26) ~af + r 2 x 2 + ... + r n x n = r , with t 1 = -x, - q 2 x 2 - ... - q n x n + q , where 

£2 Z? 

r. = n, - a i q j , i = 2, n , r = b- a x q while q i = — , r = — . Then x t = x° , 

n |_u, J 

x, = x 2 ,...,x„ = x° is a particular solution of equation (25) if and only if 
t l = = -x, - g 2 x 2 - ... - q n x° n + q , x 2 ,...,x n = x° is a particular solution of equation 

(26). 

Proof: 

x, = x[' , x 2 = x 2 ,...,x n = x°, is a particular solution of equation (25)<=> 
njXj 0 + n,x 2 + ... + a n xl = 6 <=> a { x “ + ( r 2 + a Y q 2 )x 2 + ... + (r n + a x q n )x ° = a { q + r <=> 
r 2 x 2 + ... + r n x n - n,(-x, - c/ 2 x 2 - ... - q n x n + c/) = r <=> -n^ + r 2 x 2 + ... + rx„ = r <=> 
<=>q =h°,x 2 =x 2 ,...,x n =x® is a particular solution of equation (26). 



Lemma 3. x ; =c n kj + ... + c m ] k H _ ] + r/ ( . , i = \,n , is the general solution of equation 
(25) if and only if 

(28) fj = -(c u + q 2 c 2l + ... + q„c n i )Zq - ... - (c lfl _, + q 2 c 2n _ t + ... + q„c m _ l )k n - 
~(d\ +q 2 d 2 + ... + q„ d n ) + q , 
x . = q + .. . + c jn _ x k n _, + d j , j = 2, n 
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is a general solution for equation (26). 

Proof: 

t 1 = q° = -Xj° - q 2 x j - ... - q n xll + q, x 2 = x 2 ,...,x n = x° is a particular solution of 
the equation (25) <» x, = x,° , x 2 = x 2 ,...,x n = x° is a particular solution of equation (26) 
3/cj = k° e Z,...,k n = k° eZ such that 

=c,A° +... + c iB _ 1 fe„°_ 1 +d t = x?, i = In <=> 3fej = /:[' eZ k n = k° n eZ, 

such that 

A = +... + c iB _ 1 fc„°_i +d t = x,° , i = 2,n , 

and 

6 = -( C ll + <?2 C 21 + - + WnJK ~ - " ( C l»-1 + «2 C 2»-1 + - + ~ 

(J, + q 2 d 2 + ... + q n d n ) + q = -x° - q 2 x Q 2 - ... - <?„x„ 0 + q = 6° 

Lemma 4. The linear equation 

(29) tfjXj + a 2 x 2 + ... + a n x n = 6 with |a, | = 1 has the general solution: 

x, = ~(a 2 k 2 + ... + - b)a x 

(30) « x,- = kj eZ 

i = 2, n 

Proof: 

Let’s consider x ; = x,’ , x 2 = x 2 ,...,x n = x°, a particular solution of equation (29). 
3fc 2 = x 2 , k n = x° , such that x l = (-a 2 x 2 + ... + a n x' i -b}a x = x,°, x 2 = x 2 , . . . , x n = x° . 

Lemma 5. Let’s consider the linear equation a l x l +a 2 x 2 +... + a n x n =b , with 

min k| = a x and a, = a ] q j , i = 2,n . 

a s *0 

Then, the general solution of the equation is: 
x, =-(q 2 k 2 +... + q n k n -q) 

< x ; =k t e Z 
i = 2,n 
Proof: 

Dividing the equation by a 2 the conditions of Lemma 4 are met. 

Theorem of Correctness. The preceding algorithm calculates correctly the 
general solution of the linear equation a, x, +... + a n x n = b , with not all a t = 0 . 

Proof: 

The algorithm is finite according to Lemma 1. The correctness of steps 1, 2, and 3 
is obvious. At step 4 there is always min la I as not all a, = 0 . The correctness of sub- 

a s * 0 1 1 

step 6 A) results from Lemmas 4 and 5, respectively. This algorithm represents a method 
of obtaining the general solution of the initial equation by means of the general solutions 
of the linear equation obtained after the algorithm was followed several times (according 
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to Lemmas 2 and 3); from Lemma 3, it follows that to obtain the general solution of the 
initial linear equation is equivalent to calculate the general solution of an equation at step 
6 A), equation whose general solution is given in algorithm (according to Lemmas 4 and 
5). The Theorem of correctness has been fully proven. 

Note. At step 4 of the algorithm we consider a := min I a I such that the number of 

a s * 0 1 ' 1 

iterations is as small as possible. The algorithm works if we consider a := cl | + max|n s | 

5 = 1, n 

but it takes longer. The algorithm can be introduced into a computer program. 

Application 

Calculate the integer solution of the equation: 

6x l - 12x 2 - 8x 3 + 22x 4 = 14 . 

Solution 

The previous algorithm is applied. 

1 . ( 6 , - 12 , - 8 , 22 ) = 2 

2. 2114 therefore the solution of the equation is in Z . 

3. h := 1; |2| + 1 ; dividing the equation by 2 we obtain: 

3x 4 = 6x 2 - 4x 3 + 1 lx 4 = 7 . 

4. n:=min||3|, |-6|, |-4|, |ll|} = 3, i = 1 

5. a ^ 1 

7. -6 = 3- (-2) + 0 

-4 = 3- (-2) + 2 
11 = 3-3 + 2 
7 = 3-2 + 1 



8. Xj = 2 x 2 + 2x 3 - 3x 4 + 2 — f t 

9. 



a 2 := 0 


a x := -3 


a 3 := 2 


b:= 1 


a A := 2 


x i : = h 




h := 2 


4. We have a new equation: 


-3 t l - 


0 • x, + 2x 3 + 2x 4 = 1 


a := min {-3,2,2} and 


i = 3 




5. a + 1 




7. —3 = 2- (—2) + 1 



0 = 2 -0 + 0 
2 = 2 • 1 + 0 
1 = 2-0 + 0 
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8. x 3 = 2t { + 0 ■ x 2 - x 4 + 0 - 1 2 . Substituting the value of x 3 in the value 
determined for x l we obtain: x l = 2x 2 - 5x 4 + 3/, - 2t 2 + 2 

9. flj := 1 a 3 := -2 

a 2 := 0 b := 1 
a 4 := 0 a: 3 := t 2 

h:= 3 

4. We have obtained the equation: 

1 • t 2 + 0 • x 2 - 2 ■ t 2 + 0 • x 4 = 1 , a = 1 , and i = 1 
6. (A) = -(0 • x 2 - 2t 2 + 0 • x 4 - 1) • 1 = 2t 2 + 1 

(B) Substituting the value of /, in the values of .r, and x 3 previously 
determined, we obtain: 

Xy = 2x 2 - 5.v 4 + 4 1 2 + 5 and 
x 3 = -x 4 + 3t 2 +2 

(C) x 2 := /cj , x 4 := k 2 , t 2 := k 3 , k { , k 2 , k 3 e Z 

(D) The general solution of the initial equation is: 

.Vj = 2k, - 5k 2 + Ak 3 + 5 

x 2 — k 3 

x 3 = — k 2 + 3 k 3 + 2 
x 4 = k 2 

k 4 ,k 2 , k 3 are parameters e Z 

REFERENCE 

[1] Smarandache, Florentin - Whole number solution of equations and 
systems of equations - part of the diploma thesis, University of Craiova, 
1979. 
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ANOTHER INTEGER NUMBER ALGORITHM TO SOLVE 
LINEAR EQUATIONS (USING CONGRUENCES) 



In this section is presented a new integer number algorithm for linear equation. 
This algorithm is more “rapid” than W. Sierpinski’s presented in [1] in the sense that it 
reaches the general solution after a smaller number of iterations. Its correctness will be 
thoroughly demonstrated. 



ANOTHER INTEGER NUMBER ALGORITHM TO SOLVE LINEAR 

EQUATIONS 



Let’s us consider the equation (1); (the case a n b eQ, i = 1 ,n is reduced to the 
case (1) by reducing to the same denominator and eliminating the denominators). Let 
d = If d I b then the equation does not have integer solutions, while if d\b 

the equation has integer solutions (according to a well-known theorem from the number 
theory). 

If the equation has solutions and d^ we divide the equation by d . Then d = 1 
(we do not make any restriction if we consider the maximal co-divisor positive). 

Also, 

(a) If all cl the equation is trivial; it has the general integer solution 



(b) 



x t = kf e Z, i = l,n, when b = 0 (the only case when the general solution is 
n -times undetermined) and does not have solution when b ^ 0 . 

If 3/, 1 < i < n such that cl = ±1 then the general integer solution is: 

f 



x, = -a, 



H a j k j- b 



and x „ = k r e. 



P. 




I 7=1 

V j*i 



J 



The proof of this assertion was given in [4]. All these cases are trivial, therefore 
we will leave them aside. The following algorithm can be written: 



Input 

A linear equation: 

n 

(2) '^ j a i x i = b , a n b eZ, a t ^ ±1, i = l,n , 

(=i 

with not all a,. =0 and (a 1 ,...,a n ) = l. 



Output 

The integer general solution of the equation. 

Method 

1. h := 1, p := 1 
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2. Calculate min jlrl, r = a ; (moda ), \r\< a, } and determine r and the pair 

\<i,j<n v J ' ' I J \) 

(i,j) for which this minimum can be obtained (when there are more possibilities we have 
to choose one of them). 

3. If |r| ^ go to step 4. 

If I r| = 1 , then 



x := r 



h~ Z a s X s +b 



~ a h - 



v 

f 



s = 1 



Xj := r 



ci. -r A- r - a, , 

«,-** + Z a ^+ — h 

a. J=1 a. 



(A) Substitute the values thus detennined of these unknowns in all the 
statements ( p ), p = 1 , 2 ,... (if possible). 

(B) From the last relation ( p ) obtained in the algorithm substitute in all 
relations: (p - l),(p - 2),...,(1) 

(C) Every statement, starting in order from ( p - l) should be applied the 
same procedure as in (B): then ( p - 2),.. .,(3) respectively. 

(D) Write the values of the unknowns x ; , i = l,n, from the initial 

equation (writing the corresponding integer number parameters from 
the right term of these unknowns with ), STOP. 

a — r 

4. Write the statement (p) : x j = t h — ! x i 

CI J 

5 . Assign Xj := t h h := h + 1 

a, := r p := p + 1 

The other coefficients and variables remain unchanged go back to step 2. 

The Correctness of the Algorithm 

Let us consider linear equation (2). Under these conditions, the following 
properties exist: 

Lemma 1. The set M = j |r|, r = a,(moda 7 ), 0<|r| < |a ; .|j has a minimum. 

Proof: 

Obviously M a N* and M is finite because the equation has a finite number of 
coefficients: n , and considering all the possible combinations of these, by twos, there is 
the maximum A[f (arranged with repetition) = n 2 elements. 
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Let us show, by reduction ad absurdum, that M * 0 . 

M 0 <=> a i = 0(mod af) Vz, / = 1 , n . Hence a. = 0(modrz ; ) Vz, j = l,n . Or this 
is possible only when \a\ = |a ; |, Vz, j = l,n, which is equivalent to 
(a l ,..,a n ) = a j , Vz el, n . But (a v ..,a n ) = l are a restriction from the assumption. It 
follows that \a\ = l,n, Vz e l,zz a fact which contradicts the other restrictions of the 
assumption. 

M ^ 0 and finite, it follows that M has a minimum. 

Lemma 2. If \r\ = min M , then Irl < \ci\, Vz el,zz . 

Proof: 

We assume conversely, that 3z 0 , 1 < z' 0 < zz such that r| > | | . 

Then |r| > min j \ci j |J = a f ^ 1 , 1 < j 0 < n . Let a pa , 1 < p {) < n , such that |a po | > |a ;o | and 
rz„ is not divided by ci) . 

P 0 J J 

There is a coefficient in the equation, | a jo | which is the minimum and the coefficients 

are not equal among themselves (conversely, it would mean that (a v ..,a n ) = a x = ±1 
which is against the hypothesis and, again, of the coefficients whose absolute value is 
greater that | a ijo | not all can be divided by a k (conversely, it would similarly result in 
(a l ,..,a n ) = a jo ^±1. 

We write [ ~a po / a Jo ] = q 0 e Z (integer portion), and r = a fp - q 0 a jo eZ . We have 

a po = r 0 (modrz 7o ) and 0 < |r 0 | < \a jo | < |z/ ( | < |r| . Thus, we have found an r 0 which \r 0 \ < |r| 
contradicts the definition of minimum given to |r| . 

Thus |r| < \a\ , Vz el,zz . 

Lemma 3. If |r| = minM = 1 for the pair of indices (i,j ) , then: 
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is the general integer solution of equation (2). 

Proof: 

Let x e = x°, e = l,n, be a particular integer solution of equation (2). Then 

y 

3 k s = x° eZ, se{l n}\ / } and t h = x° +— xf e Z (because o t - r = Ma j ) such 

a j 



that: 



x = r - a . 



x j= r ~ a j 



x J+ - 



a, - r 



a j j 



X a s X s +b = X i 



o , ~r o 

X , + — 5 X; 



5=1 

a , - r 



+ - 



V o , r-a,. o 

X WL + ^ = A 

5=1 ^ j 



and x s = k s = x s °, 5 e \{i,j} ■ 



Lemma 4. Let |r| ^ and (z, j ) be the pair of indices for which this minimum can 
be obtained. Again, let’s consider the system of linear equations: 



(3) 



h + rx i + X a s x s = b 



t h = x j+- 



5=1 

u, - r 



a j 



Then x e =x°, e = 1,/z is a particular integer solution for (2) if and only if x = x° e 



d — y 

e e {l,...,zz}\{j} and t h = = x] + — x j is the particular integer solution of (3). 

a j 

Proof: 

x e = x° e , e = l, n is a particular solution for (2) if and only if 



±«A =b< ^ X a s x ° s +a j 



r \ 

0 , a , -r o 
x + — x 



e=l 



5=1 



+ rx; = b <=> 



<=> af h + rx, 



+ eZ 

a. 



f + X a 5 X 5 0 = ^ and 

5=1 

se{ij} 

e e (l, n} \ { j) and t h = is a particular integer solution for (3). 



<=> X, = x„ 



Lemma 5. The previous algorithm is finite. 

Proof: 

When |r| = 1 the algorithm stops at step 3. We will discuss the case when 
r\ ^ 1 . According to the definition of r, |r| eN* . We will show that the row of 
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r - s successively obtained by following the algorithm several times is 
decreasing with cycle, and each cycle is not equal to the previous, by 1 . Let r\ be 

the first obtained by following the algorithm one time, rj ^ 1 then go to step 4, 

and then step 5. According to lemma 2, /, | < | a\, Vi = l,n. 

Now we shall follow the algorithm a second time, but this time for an 
equation in which r x (according to step 5) is substituted by a i . Again, according 

to lemma 2, the new |r| written |r 2 | will have the propriety: |r 2 | < /, . We will get 
to |r| = 1 because |r| > 1 and |r| < oo , and if |r| ^ 1 , following the algorithm once 
again we get |r| < |^ | and so on. Hence, the algorithm has a finite number of 
repetitions. 

Theorem of Correctness. The previous algorithm calculates the general 
solution of the linear equation correctly (2). 

Proof: 

According to lemma 5 the algorithm is finite. From lemma 1 it follows 
that the set M has a minimum, hence step 2 of the algorithm has meaning. When 
|r| = 1 it was shown in lemma 3 that step 3 of the algorithm calculates the general 
integer solution of the respective equation correctly the equation that appears at 
step 3). In lemma 4 it is shown that if r| ^ 1 the substitutions steps 4 and 5 

introduced in the initial equation, the general integer solution remains unchanged. 
That is, we pass from the initial equation to a linear system having the same 
general solution as the initial equation. The variable h is a counter of the newly 
introduced variables, which are used to successively decompose the system in 
systems of two linear equations. The variable p is a counter of the substitutions of 
variables (the relations, at a given moment between certain variables). 

When the initial equation was decomposed to |r| = 1 , we had to proceed in 

the reverse way, i.e. to compose its general integer solution. This reverse way is 
directed by the sub-steps 3(A), 3(B) and 3(C). The sub-step 3(D) has only an 
aesthetic role, i.e., to have the general solution under the form: x ; = ffk x ,...,k n _ 1 ) , 

i = \,n, f being linear functions with integer number of coefficients. This “if 
possible” shows that substitutions are not always possible. But when they are we 
must make all possible substitutions. 

Note 1. The previous algorithm can be easily introduced into a computer 
program. 

Note 2. The previous algorithm is more “rapid” than that of W. 
Sierpinski’s [1], i.e., the general integer solution is reached after a smaller number 
of iterations (or, at least, the same) for any linear equation (2). 

In the first place, both methods aim at obtaining the coefficient ±1 for at least 
one unknown variable. While Sierpinski started only by chance, decomposing the 
greatest coefficient in the module (writing it under the form of a sum between a 
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multiple of the following smaller coefficient (in the module) and the rest), in our 
algorithm this decomposition is not accidental but always seeks the smallest |r| 
and also choose the coefficients a, and a. for which this minimum is achieved. 

That is, we test from the beginning the shortest way to the general integer 
solution. Sierpinski does not attempt to find the shortest way; he knows that his 
method will take him to the general integer solution of the equation and is not 
interested in how long it will take. However, when an algorithm is introduced into 
a computer program it is preferable that the process time should be as short as 
possible. 

Example 1. 

Let us solve in Z 3 the equation 11 x - ly + lOz = -12 . 

We apply the former algorithm. 

1. h = l,p = l 

2. r =3, i = 3,j = 2 

3. |3| ^1 go on to step 4. 

„ 10-3 

4. (1) y = h-— — -z = t,+z 

5. Assign 

y := t, h := 2 

a, := 3 p := 2 

with the other coefficients and variables remaining unchanged, go back to 

step 2. 

2. r = -l, i = 1,7 = 3 

3. |-1| = 1 

v = — 1(— 3f 2 - (-lt x )- 12) = 3 1 2 - It j - 12 
z = —1^17 1 2 + (- It ,)■ 17 ~ (-1) + ~ 1 ~ 17 (-12)j = 11 1 2 + 42 1 , - 72 

(A) We substitute the values of x and z thus determined into the only 
statement (p) we have: 

(1) y = t x +z = — 17*2+43^-72 

(B) The substitution is not possible. 

(C) The substitution is not possible. 

(D) The general integer solution of the equation is: 

x = 3k x - 1 k 2 + 12 
< y = -Ilk, + 43L, - 72 
z = -Ilk, + 42k 2 - 72; k,,k 2 e Z 
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INTEGER NUMBER SOLUTIONS OF LINEAR SYSTEMS 



Definitions and Properties of the Integer Solution of a Linear System 

Let’s consider 

n 

(i) 'Tj a ij x j =b r 1 = l,m 

7=1 

a linear system with all coefficients being integer numbers (the case with rational 
coefficients is reduced to the same). 

Definition 1. x j = x ( ] , / = 1 , n , is a particular integer solution of (1) if 

n 

x Q j eZ, j = 1,/z and ^ a,, v) = b, i = l,m . 

7=1 

Let’s consider the functions /. : Z ,! — > Z, j = l,n , where h e N* . 

Definition 2. Xj = f j (k 1 ,...,k h ), j = 1 ,n , is the general integer solution for (1) if: 

n 

(a) Yj a ijfA k = b„ i = 1, 777 , irrespective of (k u ...,k h )e Z ; 

7 = 1 

(b) Irrespective of x . = x'j , j = \,n a particular integer solution of (1) there is 

e Zsuch that fj(k°,...,kl) = Xj,j = \,n. (In other words the general solution 
that comprises all the other solutions.) 

Property 1. 

A general solution of a linear system of m equations with n unknowns, 
r(A) = 777 < 77 , is undetermined (77 - in) -times. 

Proof: 

We assume by reduction ad absurdum that it is of order r , 1 < r < 77 - m (the case 
r = 0 , i.e., when the solution is particular, is trivial). It follows that the general solution is 
of the form: 

x \ =UnPi + - + u lr P r + v l 

X n = U nlPl + - + U nrPr + V n > U ih , V 7 G Z 
p h = parameters e Z 

We prove that the solution is undetennined ( 77 - 777 ) -times. 

The homogeneous linear system (1), resolved in r has the solution: 
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D 

x , = — ^ * 



D l 

+ J ± V 

m+1 "* *ft 



X = 



n m 

^ m+l 



d; 



^ X m+1 + "-+~ X k 



Let Xf = x°, i = l,n , be a particular solution of the linear system (1). 
Considering 

\x m+1 = D-k m+l 



x, = D ■ k„ 



we obtain the solution 



X l = D Ll-Kn + l+ ■■■+ D l'K + X l 



X m ~ D Z+\ ' K + 1 + ' • ' ■ + D 'Z ‘ 

•*>11+1 — ^ ' k/n + 1 -^m+l 



I = D-k + X , 

« n n ~ 



k . = parameters e Z 



which depends on the n - m independent parameters, for the system (1). Let the solution 
be undetennined (n- m) -times: 



x, — c,,k, + ... + C , k 

1 111 1 n—m n—m 



+ cl. 



(S 2) 



x n =c nl k l +... + u.,„ k„ _+d„ 



nn-m n-m n 



[c y ., d t e Z, kj = parameters e Z 

(There are such solutions, we have proved it before.) Let the system be: 

a, ,x, + ... + a hl x n =b x 

a m \X\ + ... + a mn x n = b rn 



X: = unknowns e Z , , b e Z . 

4 y i 

I. The case b { = 0, i = 1 ,m results in a homogenous linear system: 
a n x t + ... + a in x n = 0; z = l,m. 

(S 2 ) => a ;1 (r,f, + ... + c ln _ m k n _ m + dj)+ ... + a,.„ (<:■„*, + ... + c nn _ m k n - m + <)= 0 

0 = (a n c u + ... + a in c nl )k ] + ... + (a n c ln _ m + ... + a in c nn _ m )k n _ m + (cl 4, + ... + a in d n ) 
VC eZ 

For /:, = ... = k n _ m = 0 => a n d t +... + a m d n = 0 . 
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For k t = ... = k h _ j = k h+l = ... = k n _ m = 0 and k h = 1 => 

=> ( < Tl C ,7, + ••• + a m C 'n/ ! )+ ( a il^l + ••• + 0 => 

+ . . . + a,„c (l/l = 0, V i = 1, m, V h = l,n- m . 

The vectors 



V k = 



- ih 



h = 1, n - m 



\ C nhJ 



are the particular solutions of the system. 

V h , h = l,n-m also linearly independent because the solution is undetermined 
(n-m) -times \V 1 ,...,V n _ m } + d is a linear variety that includes the solutions of the system 



obtained from (S 2 ). 

Similarly for (Si) we deduce that 



(U u ' 



U. = 



UJ 



s = \,r 



are particular solutions of the given system and are linearly independent, because (Si) is 



undetermined (n-m) -times, and V 



is a solution of the given system. 



yvj 

Case (a) U l ,...,U r , V = linearly dependent, it follows that {f/,,...,f/ r } is a free 
sub-module of order r <n - m of solutions of the given system, then, it follows that there 
are solutions that belong tojVj ,—,V n _ m }+d and which do not belong to {[/),.. .,C/ r } , a 
fact which contradicts the assumption that (Si) is the general solution. 



Case (b) U v ...,U r , V= linearly independent. 

{C 1 ,...,C/ r }+V is a linear variety that comprises the solutions of the given system, which 
were obtained from (Si). It follows that the solution belongs to {V, ,...,V n _ m J + d and does 
not belong to {£/),... ,U r J+V , a fact which is a contradiction to the assumption that (Si) is 
the general solution. 

II. When there is an i e 1 ,m with /? ^ Othen non-homogeneous linear system 
a n x i + ... + a in x n = b v i = l,m 

(S 2 ) => a n (c n k, + ... + c ln _ m k n _ m + cl x )+... + a m (c nl k x + ... + c nn _ m k n _ m + d n )= b, 
it follows that 

= >(a,c u +... + a t „ c„, )k, +... + (a n c, + ... + a. c )k n m +(a.,cL + ... + a. d„ )= b t 

\ 11 11 in nl J 1 \ (1 1 n—m in nn—m / n—m \ /I 1 in n / 1 
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For k t = ... = k n _ m = 0 => a n d x + ... + a in d n = b x ; 

For k x = ... = k J , = £ /+l = ... = = 0 and kj = 1 

=> («,1 C I/ + - + a m c «/ ) + i a in d 1 +••• + ) = h i it follows that 

K^+-+«^.=o w . — w . 



a a d x +... + a in d n =b t 



; Vi = 1, m, Vj = l, n- m . 









-i i 



,j = l,n-m, are linearly independent because the solution (S 2 ) is 
undetermined (/r - m ) -times. 

(d x 



(? !) V-, j = \,n- m , and d = 



U J 



are linearly independent. 

We assume that they are not linearly independent. It follows that 

( s,c,, + ...+ s c, ^ 

1 11 n-m In-m 

d = s.F + ...+ s V = : 

1 1 n-m n-m 

VSiC„i + ■■■ + s n _ m c m _ m J 



Irrespective of i = 1, m : 

b x = a n d x + ... + a m d n = a n ( s x c u + ... + s n _ m c ln _ m )+... + a m (s x c nl + ... + s n _ m c nn _ m )= 

= (a n c n + ... + a in c nl )s 1 + .... + (a n c ln _ m + ... + a in c nn _ m )s n _ m = 0 . 

Then, b t = 0 , irrespective of i = 1, m , contradicts the hypothesis (that there is an i e 1, m , 
b, ^ 0). It follows that V l ,...,V n _ m ,d are linearly independent. 

{V l ,...,V n _ m } +d is a linear variety that contains the solutions of the non- 
homogeneous system, solutions obtained from (S 2 ). Similarly it follows that 
{G 1 ,...,G,}+V is a linear variety containing the solutions of the non-homogeneous 

system, obtained from (Si). 

n - m > r it follows that there are solutions of the system that belong to 



“?!” means “to prove that” 



{Vj ,...,V n _ m ] +d and which do not belong to [G x ,...,G r ]+V , this contradicts the fact that 

(Si) is the general solution. Then, it shows that the general solution depends on the n - m 
independent parameters. 

Theorem 1. The general solution of a non-homogeneous linear system is equal to 
the general solution of an associated linear system plus a particular solution of the non- 
homogeneous system. 

Proof: 
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Let’s consider the homogeneous linear solution: 

\a n x l +- + a hl x il =0 

, (AX= 0) 



a.x , + ... + ax„ = 0 

ml 1 mn n 



with the general solution: 

[x, = c u k, + ... + c. k + d, 

I 1 111 1 n—m n—m 1 



x = c ,k, + ... + c k + d 

n n 1 1 nn-m n-m n 



and 



with the general solution a particular solution of the non-homogeneous linear system 
AX = b ; 



(?!) 



x l =c u k l +... + c u _ m k n _ m +d + x i 



x =c,k*+... + c k +d +jt 

n n\ 1 nn-m n-m n n 



is a solution of the non-homogeneous linear system. 
We note: 



A = 

(vector of dimension m ), 



'a u ...a hl > 


, x = 




, b = 




, 0 = 


'0) 


\a .... a ) 

v ml mn' 




,x n J 








vOv 



1 

K = 




i i 

, c = 


(c 

C ll** 


c 'l 

* In— m 


\ i 

, d = 




\ 

,*° = 


fx^ 










'* ^ nn-m ^ 




Kd n j 




I 0 



AX = A (Ck + d + x°)= A(Ck + d)+ AX° = b + 0 = b 
We will prove that irrespective of 

*i = y° 



x = y n 

there is a particular solution of the non-homogeneous system 



27 




k x =k 1 °e Z 
k = k° eZ 

n-m n-m 

with the property: 



*i =c n k°+... + c ln kl m +d x + = y* 

X n = C ,X + - + C nn-,nK-,n + C h + X l = Yn 

,<G 

/l 

We note Y° = : 

o 

We’ll prove that those k 1 ’ e Z, j = l,n - m are those for which A(CX 0 + 0 

(there are such X° eZ because 



Xj = 0 
= 0 



is a particular solution of the homogeneous linear system and X = CK + d is a general 
solution of the non-homogeneous linear system) 



A(CK° +d + X° -Y°) = A(CK° +d) + AX° -AY 0 = 0 + b-b = 0 . 



Property 2 The general solution of the homogeneous linear system can be written 
under the form: 

(SG) 

x i = c u k l + -+c ln _ m k n _ m 

( 2 ) : 

x = c ,k, + ... + c k 

n n 11 nn-m n-m 

kj is a parameter that belongs to Z (with c/, = d 2 = ... = d n = 0 ). 

Poof: 

(SG) = general solution. It results that (SG) is undetermined (n - m ) -times. 

Let’s consider that (SG) is of the form 

X, = C, ,p, +... + C, D + d, 

1 llr 1 1 n-mr' n-m 1 

(3) : 

x = c ,p, + ...+ c p + d 

Y n n lr 1 tin- mr n-m n 

with not all d t , = 0 ; we’ll prove that it can be written under the form (2); the system has 

the trivial solution 
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x x = 0 e Z 



x= 0 eZ 



it results that there are Pj e Z, j = 1, n - m , 

k =C llA° +- + C ln- m Plm+ d l = 0 



( 4 ) 



X n = C nlPl +- + C nn-mPlm +d n =0 



Substituting p. = kj + Pj , j = \,n - m in (3) 
k. eZ) 

p-,**r p ' ez - 

Pi 

P ; e zt =>tj “ p '" pJeZ 

which means that that they do not make any restrictions. 
It results that 



*i = c nK + - + c u- m k n - m + (cn Pi + - + c u _ m pl_ m + d x 

X n = C n\K + •" + C nn-mK-m + (fnlPl + - + C m-mPl-m + d 



But 



ChiPi + - + c hn - m pLm + d h =0,h = l,n (from (4)). 
Then the general solution is of the form: 



x, — Ci , k* + . . . + c, k 

1 111 ln-m n—m 



X n ~ C nlK + ■■■+ C nn-mK-m 



k. = parameters e Z , j = l,n - m ; it results that d A = d 2 = ... = d H = 0 . 



Theorem 2. Let’s consider the homogeneous linear system: 



a ll x 1 + ... + a ln x n = 0 



a ,x + ... + a,„ x = 0 

ml 1 mn n 



r(A ) = m , (a /ll ,...,a /!n )= 1 , h = l,m and the general solution 

\x l =c u k l + ... + c l n_ m k n _ m 



X n ~ C nlK + ■••+ C nn-mK- 



then 
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( a hl ’ • • •’ a hi -\ ’ ^Ai+1 ’ • • • ’ fl /m ) ( C i 1 ’ • • • ’ C in-m ) 



irrespective of h = 1, m and i = 1, n . 

Proof: 

Let’s consider some arbitrary h e\,m and some arbitrary i el,//; 

+ - + %-!*/-! + %+!*/+! + - + = «**,- • 

Because 

( a ;jl , . . . , d hi _i , a hi+l ? • • • ? a hn ) I a hi 

it results that 

cl = (a M ,...,a hi] , a hM ,...,a hn ) x ( . 

irrespective of the value of x i in the vector of particular solutions. 

For k 2 = k i = ... = k n _ m = 0 and k x = 1 we obtain the particular solution: 



\ x i= c n =>d\c n 



For k { = k 2 = ...= k n _ m _ x = 0 and k n _ m = 1 it results the following particular 
solution: 



x = c => cl\ c ; 

i in-m in-m 5 



I "X'n ^ nn-m 



hence 



Theorem 3. 
If 



cl I c ip j = 1 ,n - m => d\(c a ,...,c in _ m ). 



x, =c r] k ] +... + c u _ m k n _ 



x = c ,k. + ... + c k 

^ n n 11 nn-m n-m 

kj = parameters e Z , c i; e Z being given, is the general solution of the homogeneous 
linear system 

a n Xj + ... + a ln x n = 0 
< : , r(A) = m < n 

a ,x, +... + a x = 0 

l ml 1 mn n 
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then (c Xj ,...,c nj ) = 1 , V/' = l,n - m . 

Proof: 

We assume, by reduction ad absurdum, that there is j 0 e\,n- m : (c ]h ,—,c njo ) = cl 

we consider the maximal co-divisor > 0 ; we reduce to the case when the maximal co- 
divisor is -d to the case when it is equal to d (non restrictive hypothesis); then the 
general solution can be written under the form: 

*1 = CnK + - + c ij 0 dk j 0 + - + Cu- m k- m 

(5) : 

x = c ,k, + ... + c . dk +... + C k 

n n\ 1 nj 0 j Q nn-m n-m 

where d = (c„. ,..,c . ), c„ =d-c ll and (c„. ,...,c . 1 = 1. 

V ‘Jo ’ ’ n Jo J ’ ‘Jo ‘Jo \ ‘Jo’ ’ “Jo ) 

We prove that 



*i = c i,„ 



*» = C “io 



is a particular solution of the homogeneous linear system. 
We’ll note: 

(k x \ 

(c n ...c, d ... c IfI m "l : 



C = \ : 



c .... c d ... c 

v nl njQ n. 



, k = k: 



x = C -k the general solution. 

a n ... a ln 

We know that AX = 0 => A(CK ) = 0 , A = : 



We assume that the principal variables are x v ...,x m (if not, we have to renumber). 
It follows that x m+l ,...,x n are the secondary variables. 

For k,=... = k. , = k , : J . 1 = ... = k , = 0 and k . =1 we obtain a particular 

1 Jo-1 J o +1 n ~ m Jo r 

solution of the system 

*1 = C ' l jo d (c ^ pi =C 0„ 

<: => 0 = A : =d-A : => A : = ()=><: 



x = c d 

I " n/o 



f cj . ^ r cj . ^ 

I WO I I Wo I 



*i =c u„ 



x„ = c„ 



is the particular solution of the system. 

We’ll prove that this particular solution cannot be obtained by 
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( 6 ) 



x = c u k, + ... + c u dk: +... + C. k =c,, 

1 11 1 lJo lo 1 n-m n-m 1 j 0 



X n = c n A + - + c nj 0 dk jo + ... + C nn _ m k n _ m = c nj0 



(V) 



X m + 1 - C m+A + - + C m+ ldkj o + ■■■ + C m+ l, n - m k n _ m ~ C, n+l jo 



x = c ,k, + ... + c ■ dk +... + C k =c 

n n\ 1 n Jo Jo nn-m n-m nj 0 



k = 

Jo 



c c c 

m+1,1 *** m+lj **' m+l,n-m 

: : 0 . : 



-h, i 



C nj 



c c d c 

^m+ 1,1 *** l " m+lj 0 *** l "'m+l,rt-m 



^,1 



: 0 . : 

C njd ••• C n,n-m 



= -gZ 
d 



(because d ^ 1). 

It is important to point out the fact that those kj = k" , j = 1 ,n - m , that satisfy the 

system (7) also satisfy the system (6), because, otherwise (6) would not satisfy the 
definition of the solution of a linear system of equations (i.e., considering the system (7) 
the hypothesis was not restrictive). From X. eZ follows that (6) is not the general 

solution of the homogeneous linear system contrary to the hypothesis); then 
(c Xj ,...,c n ^ = 1, irrespective of j = 1 ,n - m . 



Property 3. Let’s consider the linear system 

\a u x 1 + - + a hl x n = b x 

a m i x i + ••• + a mn x n = b m 
a.. ,b. eZ, r(A) = m < n , x, = unknowns e Z 
Resolved in E , we obtain 

l X l = fl( X „ + l’-’ X n) 



, x 1 ,...,x m are the main variables, 



where f are linear functions of the fonn: 



fi = 



C m + l X ,n +l +-+ C n X n+ e i 



d; 



where c‘ m+j , d t , eZ ; i = 1 ,m, j = 1 ,n - m . 
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c — 

If — eZ irrespective of i = 1 , m then the linear system has integer solution. 

d, 

Proof: 

For \<i< m, x i eZ, then fj e Z . Let’s consider 

\ X m+ 1 — U m + l^m+l 



x„ = uk 

n n n 



X \ = V m+lK+l + ••• + V nK + ~7 

c h 



m 7 . . m i , m 

X m = V m+ l k ,n + l+-+ V nK+ — 

d m 

a solution, where u m+l is the maximal co-divisor of the denominators of the fractions 



"m+J 

d : 



, i' = l, m, j = 1, n - m calculated after their complete simplification. 



C hi 

vi . = m+J " ,+J eZ is a (n - m) -times undetennined solution which depends on 



m+J 



d , 



/ 1 - 777 independent parameters (k m+v ...,k n ) but is not a general solution. 
Property 4. Under the conditions of property 3, if there is an 



i'o e 1, 777 : f in = u* +] x m+1 + ... + u’ n °x n +-f- with u'J n+j eZ, j = 1, n - m , and -f- £ Z then the 



d 



d, 



system does not have integer solution. 

Proof: 

\/x m+l ,...,x n in Z , it results that x t £ Z . 



Theorem 4. Let’s consider the linear system 

a n x 1 + ... + a ln x n =h, 

< ; 

a , jc, + ... + a x = b 

L. ml 1 mn n m 

a t] , b i eZ, Xj= unknowns eZ , r(A) = m<n. If there are indices 1 < 7i < ... < i m < n, 
i h e {1,2, .., 77 } , h = 1,777 , with the property: 
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^ 0 and 



A = 



a . ... a . 

mu mi,„ 




b x a, ... a u 



? a ...a 

m mu n 



is divided by A 



A, = 



a ... a b x 



a ... a b 

mi\ ml m-l m 



is divided by A 



then the system has integer number solutions. 
Proof: 

We use property 3 

cl = A, i = l,m; <?, = A r , h = 1 ,m 

1 l h •A'ii, 



Note 1. It is not true in the reverse case. 



Consequence 1. Any homogeneous linear system has integer number solutions 
(besides the trivial one); r(A) = m < n . 

Proof: 

l \ x = 0 : A , irrespective of h = 1, m . 

ih 

Consequence 2. If A = ±1 , it follows that the linear system has integer number 
solutions. 

Proof: 

A x : (±1) , irrespective of h = 1, m ; 

A x eZ. 

4 
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FIVE INTEGER NUMBER ALGORITHMS 
TO SOLVE LINEAR SYSTEMS 



This section further extends the results obtained in chapters 4 and 5 (from linear 
equation to linear systems). Each algorithm is thoroughly proved and then an example is 
given. 

Five integer number algorithms to solve linear systems are further given. 



Algorithm 1 . (Method of Substitution) 

(Although simple, this algorithm requires complex computations but is, 
nevertheless, easy to implement into a computer program). 

Some integer number equation are initially solved (which is usually simpler) by 
means of one of the algorithms 4 or 5. (If there is an equation of the system which does 
not have integer systems, then the integer system does not have integer systems, then 
Stop.) The general integer solution of the equation will depend on n-1 integer number 
parameters (see [5]): 



(Pi) 





i = 1 ,n , 



where all functions /.‘ 1 1 are linear and have integer number coefficients. 



This general integer number system (/;, ) is introduced into the other m — 1 
equations of the system. We obtain a new system of m - 1 equations with n - 1 unknown 
variables: 



k^, z, = l,n-l, 

which is also to be solved with integer numbers. The procedure is similar. Solving a new 
equation, we obtain its general integer solution: 

(ft) = 4 e '(c.-.eA k = mw 

where all functions f- 2t are linear, with integer number coefficients. (If, along this 

algorithm we come across an equation, which does not have integer solutions, then the 
initial system does not have integer solution. Stop.) 

In the case that all solved equations had integer system at step ( /), 1 < j <r 
( r being of the same rank as the matrix associated to the system) then: 

ip/) = f-. J) (k\ J \-Mjlj), ij = l,n-j + l, 

f (j) are linear functions and have integer number coefficients. 



Finally, after r steps, and if all solved equations had integer solutions, we obtain 
the integer solution of one equation with n-r + 1 unknown variables. 

The system will have integer solutions if and only if in this last equation has 
integer solutions. 

If it does, let its general integer solution be: 



(Pr) 





r + 1, 
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where all f. n are linear functions with integer number coefficients. 

We’ll present now the reverse procedure as follows. 

We introduce the values of k- r l> , i r = l,n - r + 1 , at step p r in the values of 





2) , K-\ = 1, n - r + 2 


from step (p r _ 1 ) . 




It follows: 




sL T 

'■K 

ii 

r? 

1 __ 
1 


+ (/«'' (kg,..., kg),-. fl'-L (kg,-, kg,)) = gg" (kg,..., kg 



l,n - r- 1 , from which it follows that g\ r 11 are linear functions with integer 
number coefficients. 

Then follows those ( p r _ 2 ) in ( p r e ) and so on, until we introduce the values 
obtained at step ( p 2 ) in those from the step (/?,). 

It will follow: 

notation g. [k x ,...,k n _ r j, i =1 ,n , with all g. most obviously, linear functions with 

integer number coefficients (the notation was made for simplicity and aesthetical 
aspects). This is, thus, the general integer solution, of the initial system. 



The correctness of Algorithm 1. 

The algorithm is finite because it has r steps on the forward way and r - 1 steps 
on the reverse, (r < +co) . Obviously, if one equation of one system does not have (integer 
number) solutions then the system does not have solutions either. 

Writing S for the initial system and 5 . the system resulted from step ( p j ) , 

l<j<r-2, it follows that passing from (/? ; ) to (p j+1 ) we pass from a system S j to a 
system S j+l equivalent from the point of view of the integer number solution, i.e. 



which is a particular integer solution of the system S J if and only if 




i j+l = l,n-j. 



is a particular integer solution of the system S j+l where 



k l=f^ 



O+l) 



(h ,—,t n _ j+l y i j+1 i, 



n~J 



Hence, their general integer solutions are also equivalent (considering these 
substitutions). Such that, in the end, resolving the initial system S is equivalent with 
solving the equation (of the system consisting of one equation) S r _ j with integer number 
coefficients. It follows that the system S has integer number solution if and only if the 
systems 5) have integer number solution, 1 < j < r - 1 . 



Example 1. By means of algorithm 1, let us calculate the integer number solution 
of the following system: 
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( 5 ) 



r . 5x - ly - 2z + 6w = 6 
[-4a: + 6y - 3z + 1 lw = 0 

Solution: We solve the first integer number equation. We obtain the general 
solution (see [4] or [5]): 

x = t x + 2 t 2 

y = h 

z = -t 1 + 5 t 2 + 3t 3 - 3 
w = t 3 

where t v t 2 , r 3 eZ. 

Substituting in the second, we obtain the system: 

(SJ 5t l -23t 2 +2t 3 + 9 = 0. 

Solving this integer equation we obtain its general integer solution: 

h = K 

( p 2 ) < t 2 = k l + 2 k 2 + 1 

t 3 = 9/r| + 23 k 2 + 7 

where k v k 2 eZ. 

The reverse way. Substituting ( p 2 ) in (p l ) we obtain: 

x = 3Atj + 4 k 2 + 2 
y = k i 

z = 3 1 Arj + 19k 0 + 23 
w = 9k\ + 23 k 2 + 7 

where k v k 2 eZ, which is the general integer solution of the initial system (5) . Stop. 

Algorithm 2. 

Input 

A linear system (1) without all a tj = 0 . 

Output 

We decide on the possibility of an integer solution of this system. If it is possible, 
we obtain its general integer solution. 

Method 

1 . t = 1 , h = 1 , p = 1 

2. (A) Divide each equation by the largest co-divisor of the coefficients of the 
unknown variables. If you do not obtain an integer quotient for at least one equation, then 
the system does not have integer solutions. Stop. 

(B) If there is an inequality in the system, then the system does not have integer 
solutions. Stop. 

(C) If repetition of more equations occurs, keep one and if an equation is an 
identity, remove it from the system. 
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3. If there is (L, z n ) such that \a. =1 then obtain the value of the variable 

v u J u 7 I l oJo | 

x jo from equation z 0 ; statement (T t ) . 

Substitute this statement (where possible) in the other equations of the system and 
in the statement (7V, ) , (H h ) and (P p ) for all z, h , and p . Consider /:=/ + !, remove 



equation (z 0 ) from the system. If there is no such a pair, go to step 5. 

4. Does the system (left) have at least one unknown variable? If it does, consider 
the new data and go on to step 2. If it does not, write the general integer solution of the 
system substituting k v k 2 ,... for all variables from the right term of each expression 
which gives the value of the unknowns of the initial system. Stop. 

5. Calculate 



a = mm \\r\a 

UJ1J2 



= r (moda i:; j, 0< 



r < 



a. 



'h 



and determine the indices z, /, , j 2 as well as the r for which this minimum can be 
calculated. (If there are more variables, choose one arbitrarily.) 

a tj - r 

6. Write: x j =t h — 1 x n , statement (H h ) . Substitute this statement (where 

72 a u 11 

l Jl 

possible in all the equations of the system and in the statements (T I ), (H lt ) and ( P p ) for 



all t, h , and p . 

7. (A) If a ^ 1 , consider x h := t h , h := h + 1 , and go on to step 2. 

(B) If a = 1 , then obtain the value of Xj from the equation (i) ; statement 
(P p ) . Substitute this statement (where possible in the other equations of the 
system and in the relations (7) ), (H lt ) and ( P p , j for all t, h , and p . 



Remove the equation (z) from the system. 

Consider h := h + 1, p := p + 1 , and go back to step 4. 



The correctness of algorithm 2. Let consider system (1). 



Lemma 1. We consider the algorithm at step 5. Also, let 

M ={H> a ij \ = r(moda ij2 ), 0 < |r| < a ih , i, j v j 2 =1,2,3,... 



Then M ^ 0 . 

Proof: 

Obviously, M is finite and M c fL . Then, M has a minimum if and only if 
M^0. We suppose, conversely, that M = 0 . Then 
% = 0(mod a lh ), V i,j 1 ,j 2 . 



It follows as well that 



That is 



a ij 2 = 0(mod a ih ), V z, , j 2 . 



\ a ij\ = \ a ij 2 \’ v iJvji- 

We consider an z 0 arbitrary but fixed. It is clear that 
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K,p •••’%) : 

(because the algorithm has passed through the sub-steps 2(B) and 2(C). But, because it 
has also passed through step 3, it follows that 

kb 1 * v/* 

but as it previously passed through step 2(A), it would result that 

kb 1 ’ v/- 

This contradiction shows that the assumption is false. 



Lemma 2. Let’s consider a . . = r( mod a,, ) . Substitute 

IqJ\ v lh ' 



IJ2 ' 



a i i ~ r 

x h =t x h 



a 



l oJi 



in system (A) obtaining system (B). Then 

Xj = x°, j = 1 ,n 

is the particular integer solution of (A) if and only if 



n n a i i ~ Y 

A = A . J * h and h = X U 



a 



l oJi 



is the particular integer solution of (B). 



Lemma 3. Let a, ^ and obtained at step 5. 

Then 0 < a 2 < a, 

Proof: 

It is sufficient to show that a, < |a,J, V i, j because in order to get a 2 , step 6 is 
obligatory, when the coefficients if the new system are calculated, a, being equal to a 
coefficient form the new system (equality of modules), the coefficient on (i 0 /, ) . 



Let a t j with the property a ln /n < n, . 



Hence, > 



l oJ 0 



a. ■ 

l oJ 



= mm 



a. 



IqJ 



. Let a, . with a ; . > a i; : there is such an element 

*0 Jr lr> ' ■ 11 - 



*OJs 



because a, , is the minimum of the coefficients in the module and not all a , , j = l,n 

^OJr” m / 7 •/ 



| 



are equal (conversely, it would result that (a. . ) ~ a. ., V/el,r, the algorithm 

passed through sub-step 2(A) has simplified each equation by the maximal co-divisor of 
its coefficients; hence, it would follow that <a, = f Y/' = l,n , which, again, cannot be 
real because the algorithm also passed through step 3). Out of the coefficients a t . . we 
choose one with the property a, ^ Ma : there is such an element (contrary, it would 

*0 J Sq J m 



result (a. .,...,n.J~ 




but the algorithm has also passed through step 2(A) and it 



would mean that | a ioj | = 1 which contradicts step 3 through which the algorithm has also 
passed). 
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Considering q 0 = 



ci / a 

l 0JSa l 0Jm 



e Z and r = a. - q n a. e Z , we have 

l oJs 0 1U l n./ m 



l oJs„ 



= r 0 (modn ( . . ) and 0 < \r 0 \ < cl . < cl . < a, . We have, thus, obtained an r 0 with 



r 0 < a, , which is in contradiction with the very definition of a x . Thus a x < a lj , Vi, j . 



Lemma 4. Algorithm 2 is finite. 

Proof: 

The functioning of the algorithm is meant to transform a linear system of m 
equations and n unknowns into one of m, x n, with m, < m , «, < n , thus, successively 
into a final linear equation with n - r + 1 unknowns (where r is the rank of the 
associated matrix). This equation is solved by means of the same algorithm (which works 
as [5]). The general integer solution of the system will depend on the n - 1 integer 
number independent parameters (see [6] - similar properties can be established also the 
general integer solution of the linear system). The reduction of equations occurs at steps 
2, 3 and sub-step 7(B). Step 2 and 3 are obvious and, hence, trivial; they can reduce the 
equation of the system (or even put an end to it) but only under particular conditions. The 
most important case finds its solution at step 7(B), which always reduces one equation of 
the system. As the number of equations is finite we come to solve a single integer number 
equation. We also have to show that the transfer from one system m l x il to another 
m i+] x n i+l is made in a finite interval of time: by steps 5 and 6 permanent substitution of 
variables are made until we to n = l(we to a = 1 because, according to lemma 3, all 
a - s are positive integer numbers and form a strictly decreasing row). 

Theorem of correctness. 

Algorithm 2 correctly calculates the general integer solution of the linear system. 

Proof: 

Algorithm 2 is finite according to lemma 4. Steps 2 and 3 are obvious (see also 
[4], [5]). Their part is to simplify the calculations as much as possible. Step 4 tests the 
finality of the algorithm; the substitution with the parameters k v k 2 ,... has 
systematization and aesthetic reasons. The variables t, h,p are counter variables (started 
at step 1) and they are meant to count the statement of the type T,H,P (numbering 
required by the substitutions at steps 3, 6 and sub-step 7(B); h also counts the new 
(auxiliary) variables introduced in the moment of decomposition of the system. The 
substitution from step 6 does not affect the general integer solution of the system (it 
follows from lemma 2). Lemma 1 shows that at step 5 there is always a , because 
0 ^ M c N* . 

The algorithm performs the transformation of a system m l x n i into another 
m i+l x n i+l , equivalent to it, preserving the general solution (taking into account, 
however, the substitutions) (see also lemma 2). 

Example 2. Calculate the integer solution of: 
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\2x -ly + 9z = 12 
- 5y + Sz + 10 w = 0 
0z + 0 w = 0 
15.x: + 21z + 69w = 3 

Solution: 

We apply algorithm 2 (we purposely selected an example to be passed through all 
the steps of this algorithm): 

1 . t = 1, h = 1, p = 1 

2. (A) The fourth equation becomes 5x + lz + 23w = 1 

(B) - 

(C) Equation 3 is removed. 

3. No; go on to step5. 

5. a = 2 and i = 1, j l = 2, j 2 = 3 , and r = 2 . 

6. z = t l + y, the statement (//,). Substituting it in the 

\2x-2y + 9t x =12 
< 3y + 9t l + 10 w = 0 

5x +7y + 7t 1 +23w =1 

7. a ^ 1 consider z = t v h := 2 , and go back to step 2. 

2.- 

3. No. Step 5. 

5. a = 1 and i = 2, j t = 4, j 2 = 2 , and r = 1 . 

6. y = t 2 - 3w , the statement (H 2 ) . Substituting in the system: 

-12x + 2t 2 +9tj -6w = 12 
* 3?2 + + w = 0 

5.x: +7 1 2 +7fj +2w = l 

Substituting it in statement (//,), we obtain: 
z. = t { +t 2 ~ 3vv , statement (H ] ) . 

7. w = -3 1 2 - 8r, statement ( 7) ) . 

Substituting it in the system, we obtain: 

f— 12.x: — 207 2 +57fj = 12 
[ 5.x: + 1 2 - 9t x =1 

Substituting it in the other statements, we obtain: 
z = 10t 2 +25q, (/fj" 

y = 10t 2 +24q, (// 2 )" 
h := 3, p := 2 , and go back to step 4. 

4. Yes. 

2.- 

3. t 2 = 1 - 5.x: + 9^ , statement (7) ) . 

Substituting it (where possible) we obtain: 
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{-1 12.x; + 237r, = -8 (the new system); 
z= 10-50^ + 115^, (H x ) m 
y = 10-50x + 114?!, {H 2 ) n 
w = -3 + 15.x; + 35tj, (P x )' 

Consider ? := 2 go on to step 4. 

4. Yes. Go back to step 2. (From now on algorithm 2 works similarly with that 
from [5], with the only difference that the substitution must also be made in the 
statements obtained up to this point). 

2.- 

3. No. Go on to step 5. 

5. a = 13 (one three) and i = 1, j 1 = 2, j 2 = 1 , and r = 13 . 

6. x = t 3 + 2t x , statement (H 3 ). 

After substituting we obtain: 

{-1 12r 3 + 13h = -8 (the system) 

2 = 10-50/3+15?!, (H^; 
y = 10 — 5 0G +14?!, ( H 2 ) 
w = —3 + 1 5? 3 -5h, (/>!)"; 

? 2 = 1 — 5? 3 — h (7))'; 

7. a; := ? 3 , h := 4 and go on to step 2. 

2.- 

3. No, go on to step 5. 

5. a = 5 and i = 1, y, = 1, j 2 = 2 and r = 5 

6. ?j = ? 4 + 9? 3 , statement (// 4 ) . 

Substituting it, we obtain : 

5?3 + 13? 4 = -8 (the system). 

z= 10 + 85? 3 +15? 4 , (H x ) v - 
y = 10 + 7 6? 3 +14? 4 , ( H 2 ) ,v ; 
x= 19? 3 + 2? 4 , (77 3 )'; 
w = -3-30? 3 - 5? 4 , (/))'"; 

? 2 = 1 — 14? 3 - ? 4 (7J)"; 

7. ?j := ? 4 ;/r := 5 and go back to step 2. 

2.- 

3. No. Step 5. 

5. a = 2 and i = 1, j x = 2, j 2 = 1 and r = -2 . 

6. ? 3 = ? 5 - 3? 4 statement (// 5 ) . After substituting, we obtain: 

5? 5 — 2? 4 = -8 (the system). 
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z = 10 + 85t 5 -240t 4 , 
y = 10 + 76t 5 -214t 4 , (// 2 ) v ; 
x = 19t 5 -55t 4 , {H,) iv - 

w = -3-30t 5 + 85t 4 , (^) /y ; 

C = —1 — 1 4r 5 +4\t 4 , (Tj)'"; 

f, = 9? 5 +26? 4 , 

7. t 3 := t 6 ,/r := 6 and go back to step 2. 

2.- 

3. No. Step 5. 

5. a = 1 and i = 1, j l = 2, j 2 , r = 1 . 

6. t 4 = t 6 + 2t s statement (H 6 ) . After substituting, we obtain: 

r 5 — 2r 6 = —8 (the system) 
z = 10 — 395f 5 -240f 6 , {H x ) vn - 
y = 10-392t 5 -214f 6 , ( H 2 ) ,v ; 
x = —9lt s —55t 6 , (H 3 ) ; 

w = -3 + 1406, + 85t 6 , (i>) v ; 

t 2 = l + 68t 5 +4h 6 , (T l ) IV ; 

h= -43t 5 -26t 6 , (H 4 )"; 

t 3 — — 5t 5 —3t 6 , (7/j), 

7. t 5 = 2r 6 - 8 statement (P 2 ). Substituting it in the system we obtain: 0=0. 
Substituting it in the other statements, it follows: 

z = -1030/g + 3170 
y = -918f 6 + 2826 
x = -237 t 6 + 728 
w= 365t 6 -1123 
t 2 = 177t 6 -543 
f, =112t 6 +344 

> statements of no importance. 

? 3 = 13t 6 + 40 
t 4 = 5 1 6 - 16 

Consider h := l,p := 3 , and go back to step 4. t 6 eZ 

4. No. The general integer solution of the system is: 
x = -237k, +728 

y = -918k, +2826 
z = 1030k; +3170 
w= 365k, -1123 
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where k x is an integer number parameter. 

Stop. 

Algorithm 3. 

Input 

A linear system (1) 

Output 

We decide on the possibility of an integer solution of this system. If it is possible, 
we obtain its general integer solution. 

Method 



1 . Solve the system in R" . If it does not have solutions in R" , it does not have 
solutions in Z” either. Stop. 

2. f = l,t = l,h = l,g = l 

3. Write the value of each main variable x i under the form: 

( E f. i) : A = ZV; + + [Z r a x j + f j 7 A < 

with all q tj , q t , r tj , r n A ; in Z such that all |^.|<|A / |, A,. ^0, |/;. | < |A f . | (where all x j of 

the right term are integer number variables: either of the secondary variables of the 
system or other new variables introduced with the algorithm). For all i , we write 



4. (E f ^ TyXj - r ijf Y f ' , + r = 0 where (Y f ) are auxiliary integer number 

j 

variables. We remove all the equations (F f which are identities. 

5. Does at least one equation \F f ^ exist? If it does not, write the general 



integer solution of the system substituting k x ,k 2 ,... for all the variables from the right 
tenn of each expression representing the value of the initial unknowns of the system. 
Stop. 

6. (A) Divide each equation (F f by the maximal co-divisor of the coefficients 
of their unknowns. If the quotient is not an integer number for at least one i 0 the system 
does not have integer solutions. Stop. 

(B) Simplify -as in m - all the fractions from the statements (E f ; . 

7. Does r. () /n exist having the absolute value 1? If it does not, go on to step 8. If it 



does, find the value of x t from the equation \F f io write (T t ) for this statement, and 
substitute it (where it is possible) in the statements (£/„) (r-) (h„> (o,) for all 

i, t, h and g . Remove the equation (F f h Consider / := / + 1, t := t + 1 , and go back 

to step3. 

8. Calculate 
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a = min \ \r \ , r jf = r(mod r u ), 0 < \r\ < 

‘J1J2 



IJl V i/2 ' 

and determine the indices i m , j 2 as well as the r for which this minimum can be 

obtained. (When there are more variables, choose only one). 

a j . -r , 

9. (A) Write x, = z h — x j , where z h is a new integer variable; statement 

J2 a. . 

JmJl 

(«.)• 

(B) Substitute the letter (where possible) in the statements 
(£, ,) (F/,;) (T t ), (H h ,), (G g ) for all i, t, h and g. 

(C) Consider h := h + 1 . 

10. (A) If a 1 go back to step 4. 

(B) If a = 1 calculate the value of the variable v ; from the equation (F f ( 
relation (g‘ Substitute it (where possible) in the statements (£,,) (T r ), (H„), (G g „) 
for all i, t, h , and g . Remove the equation (F f J. Consider g := g + 1, f := f + 1 and go 
back to step 3. 

The correctness of algorithm 3 



Lemma 5. Let i be fixed. Then 



Z c-c • r 



\j= n i 



A ( . (with all r, r n A ; , n v n 2 being 



integers, n, < n 2 , A, ^ 0 and all Xj being integer variables) can have integer values if 
and only if r ilh , A,.) | r r 

Proof: 

The fraction from the lemma can have integer values if and only if there is a 
z e Z such that 

f n-, h n 2 

\^i=Z^Y j r jj x j -A l z + r i =0, 






\J= n 1 



J=n 1 



which is a linear equation. This equation has integer solution <=> (f,,, , A, ) f • 



Lemma 6. The algorithm is finite. It is true. The algorithm can stop at steps 1,5 or 
sub-steps 6(A). (It rarely stops at step 1). 

One equation after another are gradually eliminated at step 4 and especially 7 and 
10(B) (F f ) - the number of equation is finite. 

If at steps 4 and 7 the elimination of equations may occur only in special cases 
elimination of equations at sub step 10 (B) is always true because, through steps 8 and 
9 we get to a = 1 (see [5]) or even lemma 4 (from the correctness of algorithm 2). 
Hence, the algorithm is finite. 

Theorem of Correctness. 

The algorithm 3 correctly calculates the general integer solution of the system (1). 
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Proof: 

The algorithm if finite according to lemma 6. It is obvious that the system does 
not have solution in R" it does not have in Z" either, because Z" cz R' ! (step 1). 

The variables /, t, h, g are counter variables and are meant to number the 
statements of the type E, F, T , H and G , respectively. They are used to distinguish 
between the statements and make the necessary substitutions (step 2). 

Step 3 is obvious. All coefficients of the unknowns being integers, each main 
variable x t will be written: 



f \ 



x, = 



z 



qjXj+Ci 



I A.. 



V s J 

which can assume the form and conditions required in this step. 

Step 4 is obtained from 3 by writing each fraction equal to an integer variable 
Y f i (this being x ; eZ). 



Step 5 is very close to the end. If there is no fraction among all [E f A it means 



that all main variables x, already have values in Z , while the secondary variables of 
the system can be arbitrary in Z , or can be obtained from the statements T, H or G 
(but these have only integer expressions because of their definition and because only 
integer substitutions are made). The second assertion of this step is meant to 
systematize the parameters and renumber; it could be left out but aesthetic reasons 
dictate its presence. According to lemma 5 the step 6(A) is correct. (If a fraction 
depending on certain parameters (integer variables) cannot have values in Z , then the 
main variable which has in the value of its expression such a fraction cannot have 
values in Z either; hence, the system does not have integer system). This 6(A) also 
has a simplifying role. The correctness of step 7, trivial as it is, also results from [4] 
and the steps 8-10 from [5] or even from algorithm 2 (lemma 4). 

Ther initial system is equivalent to the “system” from step 3 (in fact, (E f ) as 

well, can be considered a system) Therefore, the general integer solution is preserved 
(the changes of variables do not prejudice it (see [4], [5], and also lemma 2 from the 
correctness of algorithm 2)). From a system m i x n t we form an equivalent system 
m i+1 x n j , with m j+1 < m i and n i+l < n i . This algorithm works similarly to algorithm 
2 . 



Example 3. Employing algorithm 3, find an integer solution of the following 
system: 

3x, + 4x 2 + 22x 4 - 8x 5 = 25 
< 6x l + 46x 4 - 12 x 5 = 2 

4x 2 + 3x 3 - x 4 + 9x s = 26 

Solution 

1. Common resolving in R 3 it follows: 
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*1 = 



*2 = 



23x 4 - 6x 5 - 1 
-3 

x 4 + 2x 5 + 24 



llx 5 +2 



2. f = l,t = l,h = l,g = l 



3. 



-7 x 4 + 2x 5 4 

-3 

6 + fl±3x L 



x 3 = 



. X 5 + 2 

- 4jc, + — 



(£.,.) 

( £ ..0 

(£.3) 



4. 



2.^4 + 3 v n -1 = 0 ( F u ) 

-*•4 4- 2x 5 — 4 y 12 —0 ( F 1i2 ) 

x 5 -3^13 4-2 = 0 (F 1>3 ) 

5. Yes. 

6 . - 

7. Yes: |r 35 1 = 1 . Then x 5 = 3>’ l3 -2 the statement {T ] ) . Substituting it in the 
others, we obtain: 



*2 = 



x 3 = 



■ 4 + 2x * ~ 1 
-3 


(£...) 


6| h+ 6 j 13 - 4 

4 


(£,. 3 ' 


+ 8+ 3 Ti3-2 + 2 

3 


(£,.3 



Remove the equation (Tj 3 ). 

Consider / := 2, t := 2 ; go back to step 3. 

2x. - 1 



Xj = -7 x 4 4- 6j 13 - 4 4- • 



-3 



x 3 = 



j, i3+5 + ii±?2k 

- 1 ly l 3 4- 8 



k.) 

C. 3) 
C. 3) 
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4 . 



2*4 +3v 21 -1 = 0 (F 2>1 ) 

* 4 + 2v 13 ~4y 22 —0 (^2,2) 

5. Yes. 

6 . - 

7. Yes |r 24 | = 1 . We obtain x 4 = -2 v l3 + 4y 22 statement (T 2 ) . Substituting it in 
the others we obtain: 

. r ,=- 28 4 , 2 + 20 ^„ ( £j ] )' 

* *2 = 3^22 + 3Y+5 (£ 2 , 2) 

*3 ~ -11*13+8 (£2,3)' 

Remove the equation (f 2 2 ) 

Consider / := 3, t := 3 and go back to step 3. 

3. 

*1 = ~22y 13 + 30v 22 + 2:Vl3 + ^ 22 1 (£3^) 

" x 2 — kl3 + Y 22 + 5 (£2, 2 ) 

x 3 = — 11^13+8 (^3,3) 

4. 2v 13 +2v 22 + 3y 31 -1 = 0 (C. ,) 

5. Yes. 

6 . - 

7. No. 

8. a = 1 and i m = 1, j 1 = 31, j 2 = 22 , and r = 1 . 

9. (A) _y 22 = z, - y 3l (statement (//, )). 

(B) Substituting it in the others we obtain: 

*1 =~ 22yi3-30z 1 +30y 31 -4+ 2};i3+2 ^ 3 2y$1 1 ( £ 3 , , ) ’ 

<a 2 = y 13 + Zj -y 31 +5 (E 32 )' 

a 3 =— llyi3 +8 (^3,3) 

2y 13 + 2zj + y 31 - 1 = 0 (£ 3 j)’ 

*4 =-2yi 3 +4z 1 -4y 13 (T 2 )' 

(C) Consider h := 2 

10. (B) y l3 = 1 - 2y ]3 - 2z, , statement (G, ). 

Substituting it in the others we obtain: 
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x t = -40 y, 3 -92z,+27 


(+,)" 




x 2 = 3y 13 + 3zj + 4 


ted" 




x 3 = -1 ly l3 +8 


ted" 




x 4 =6^3+12^ -4 


(T,Y 




3^22 = 2 Al3 + 3 ^l - 1 


ted' 




Remove equation (f 3 , ). 






Consider g := 2, f := 4 and go back to step 3. 




3. 


x l = -40y u -92z l +27 


(£.4) 




< 


*2 = 3y 13 + 3zj +4 


ted 






x 3 =-lly 13 +8 


te.d 




4. - 








5. No. The general solution of the initial system is: 




= -40k, -92 k 2 + 27, 


from ( 


£<d 




x 2 = 3k, + 3 k 2 + 4, 


from | 


ted 


< 


ll 

l 

+ 

00 


from 


(£ 43 ) 




x 4 = 6kj + 12C -4, 


from 


+2)' 




x 5 = 3 kj -2, 


from 


h) 



where k v k 2 eZ. 



Algorithm 4 
Input 

A linear system (1) with not all a tj = 0 . 

Output 

We decide on the possibility of integer solution of this system. If it is possible, we 
obtain its general integer solution. 

Method 

1. h = l, v = l. 

2. (A) Divide every equation i by the largest co-divisor of the coefficients of the 
unknowns. If the quotient is not an integer for at least one i 0 then the system 
does not have integer solutions. Stop. 

(B) If there is an inequality in the system, then it does not have integer 
solutions 

(C) In case of repetition, retain only one equation of that kind. 
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(D) Remove all the equations which are identities. 

3. Calculate a = m in 1 , a tj ■*- o| and determine the indices i 0 ,j 0 for which this 

minimum can be obtained. (If there are more variables, choose one, at 
random.) 

4. If a * 1 go on to step 6. 

If a = 1 , then: 

(A) Calculate the value of the variable Xj from the equation i 0 note this 
statement (V v ). 

(B) Substitute this statement (where possible) in all the equations of the 
system as well as in the statements (V v (H h ), for all v and h . 

(C) Remove the equation i 0 from the system. 



(D) Consider v := v + 1 . 

5. Does at least one equation exist in the system? 

(A) If it does not, write the general integer solution of the system substituting 
k v k 2 ,... for all the variables from the right term of each expression 
representing the value of the initial unknowns of the system. 

(B) If it does, considering the new data, go back to step 2. 

6. Write all a ; ^ /' and b : under the form : 

l oJ J J u *0 



a i , = a £ , cr . + r . , with r . <\a 

l o J l oJo ll o J l oJ l 0 J 



b: , = a,. , q i + r , with 

l oJ 1 qJq ± l 0 '0 



l oJ | 



*0 Jo 



7. Write x jo = - £ q kj x j + q, n + t h , statement (H h ). 

j^Jo 

Substitute (where possible) this statement in all the equations of the system as 
well as in the statement (V v ). (//,,), for all v and h . 

8. Consider 

x jo := t h , h :=h + 1, 



a := r 

'o J >oJ 



a. 



l oJo l oJo 

and go back to step 2 



J ^ J o> 

:= - a L l > b. := +r. , 



The correctness of Algorithm 4 

This algorithm extends the algorithm from [4] (integer solutions of equations to 
integer solutions of linear systems). The algorithm was thoroughly proved in our previous 
article; the present one introduces a new cycle - having as cycling variable the number of 
equations of system - the rest remaining unchanged, hence, the correctness of algorithm 
4 is obvious. 
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Discussion 

1 . The counter variables h and v count the statements H and V , respectively, 
differentiating them (to enable the substitutions); 

2. Step 2 ((A)+(B) + (C)) is trivial and is meant to simplify the calculations (as 
algorithm 2); 

3. Sub-step 5 (A) has aesthetic function (as all the algorithms described). 
Everything else has been proved in the previous chapters (see [4], [5], and 
algorithm 2). 

Example 4. Let us use algorithm 4 to calculate the integer solution of the 
following linear system: 

[ 3xj - 7x 3 + 6x 4 = -2 

[4xj + 3 x 2 + 6x 4 -5x 5 =19 

Solution 

1. h = 1, v = l 

2 . - 

3. a = 3 and i = 1, j = 1 

4. 3 + 1 . Go on to step 6. 

6. Then, 

-7 = 3- (-3) + 2 
6 = 3-2 + 0 
-2 = 3 • 0 - 2 

7. Xj =3x 3 -2x 4 +h statement (//,). Substituting it in the second equation we 
obtain: 

4fj + 3 x 2 + 12x 3 - x 4 - 5 x 5 = 19 

8. x, := f, , h := 2, a l2 := 0, a 13 := +2, a 14 := 0, a n := +3, b := -2 . 

Go back to step 2. 

2. The equivalent system was written: 

f 3 fj + 3 x 3 = -2 

[4q + 3 x 2 + 12x 3 - x 4 - 5x 5 =19 

3. a = 1, i = 2, j = 4 

4. 1=1 

(A) Then: x 4 = 4/, +3x 2 +12x 3 -5x 5 -19 statement (k, ). 

(B) Substituting it in (//,), we obtain: 

Xj = -7q - 6 x, - 21x 3 + 10x 5 + 38, (//, ) 

(C) Remove the second equation of the system. 

(D) Consider: v := 2 . 

5. Yes. Go back to step 2. 

2. The equation +3q +2x 3 = -2 is left. 

3. a = 2 and i = 1, j = 3 
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4. 2 ^ 2 , go to step 6. 

6. 

+3 = +2 • 2 - 1 
-2 = +2(-l) + 0 

7. x 3 = -2/, + 1 2 - 1 statement (/7 2 ). 

Substituting it in ( H l )' , (V, ), we obtain: 

Xj = 35q -6 x 2 -21t 2 +10x; +59 ( H 
x 4 = -20q + 3x, + 12t 2 -5x 5 -31 (V, )' 

8. x 3 := t 2 , h := 3, a n :=- 1, a 13 := +2, ^1=0 , (the others being all = 0). Go back 
to step 2. 

2. The equation -5f, + 2t 2 = 0 was obtained. 

3 . a = 1 , and i = 1, j = 1 

4. 1=1 

(A) Then t l = 2t 2 statement (V 2 ). 

(B) After substitution, we obtain: 
jcj =49t 2 -6.r 2 + 10.r 5 +59 

x 4 = -28 1 2 + 3x 2 - 5x s - 3 1 (Vj ) " ; 

x 3 =-3t 2 (H 2 )'; 

(C) Remove the first equation from the system. 

(D) v := 3 

5. No. The general integer solution of the initial system is: 

= 49 kj - 6 k 2 + 10k 3 + 59 

x 2 = k 2 

x 3 = -3kj - 1 

x 4 = -28 k 4 + 3 k 2 - 5k 3 - 3 1 
x 5 = k 3 

where (k 4 , k 2 , k 3 ) e Z 3 . 

Stop. 

Algorithm 5 
Input 

A linear system (1) 

Output 

We decide on the possibility of an integer solution of this system. If it is 
possible, we obtain its general integer solution. 

Method 

1. We solve the common system in M" , then it does not have solutions in R" , 
then it does not have solutions in Z" either. Stop. 

2. f = l, v = l, h = \ 

3. Write the value of each main variable x i under the form: 
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( E f.i\ :x i =Z«-^' + X'k V / ■ r > 



V j 



/ A 



with all q Lj , q n r tj , r, A,, from Z such that all ^ <|A ( .|, r < | A . | , A ( . ^ 



4 . 



5. 



6 . 



7. 

8 . 



9. 



(where all x ] -S of the right term are integer variables: either from the 

secondary variables of the system or the new variables introduced with the 
algorithm). For all i , we write r ; . = A ( 

(e f ^ r rj x j - r jJf y fi + r = 0 where (j/,,) are auxiliary integer variables. 

j 

Remove all the equations \F f t j which are identities. 

Does at least one equation (F f ( ) exist? If it does not, write the general integer 

solution of the system substituting k v k 2 ,... for all the variables of the right 
number of each expression representing the value of the initial unknowns of 
the system. Stop. 

(A) Divide each equation (F f ( ) by the largest co-divisor of the coefficients of 

their unknowns. If the quotient is an integer for at least one i 0 then the system 
does not have integer solutions. Stop. 

(B) Simplify - as previously ((A)) all the functions in the relations ( E f i ^ . 

Calculate a = min j|r|, r ^ o| , and determine the indices i 0 , j 0 for which this 

minimum is obtained. 

If a ^ 1 , go on to step 9. 

If a = 1 , then: 

(A) Calculate the value of the variable x Jo from the equation (F f ( ^ write 
(V v ) for this statement. 

(B) Substitute this statement (where possible) in the statement (E f i j, 
(V v+1 ), (H h ), for all i, v , and h . 

(C) Remove the equation (E f /). 

(D) Consider v := v + 1, / := / + 1 and go back to step 3. 

Write all r j * /' and r. under the fonn: 

hJ J l 0 

r. . = A . • q. . + r . , with I r. . < I A. I ; 

l 0 J *0 ll oJ l 0 J I l oJ | I l \ 

r. . = A. • q. + r. , with I r I < I A. I . 

l oJ l o ll o l o | | I 1 1 



10. (A) Write x Jo = - £ q^x] + q k> + t h statement ( H h ) . 

J^Jo 



(B) 



Substitute this statement (where possible) in all the statements ( E f 

fy,) (K). (».-,)• 
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(C) Consider h := h + 1 and go back to step 4. 

The correctness of the algorithm is obvious. It consists of the first part of 
algorithm 3 and the end part of algorithm 4. Then, steps 1-6 and their correctness were 
discussed in the case of algorithm 3. The situation is similar with steps 7-10. (After 
calculating the real solution in order to calculate the integer solution, we resorted to the 
procedure from 5 and algorithm 5 was obtained). This means that all these insertions 
were proven previously. 



Example 5 

Using algorithm 5, let us obtain the general integer solution of the system: 
f 3x l + 6x 3 + 2x 4 = 0 

| 4x 2 - 2x 3 - lx 5 = -1 



Solution 

1. Solving in M 5 we obtain: 



A 






-6.r 3 - 2x 4 
3 

-2x 3 + lx 5 - 1 

4 



2 . 

3. 



4. 



5. 

6 . 

7. 

8 . 

9. 



/ = 1, v = 1, h = 1 

(e ii ):x 1 = 2x 3 + ~y l 
s \ 2x 3 + 3.w 

\E\,2 )■ X 2 = X 5 Z 



(F u ):-2.x 4 -3y u =0 

(^1,2 ) : 2x 3 + 3-^5 _ 4 A 2 _ 1 = 0 
Yes 



i = 2 and i 0 = 2, j 0 = 3 
2^1 

3=21+1 
-4 = 2- (-2) 



1 



-1 = 20-1 

10. x 3 = -x, + 2y 12 + tj statement (//, ). After substitution: 



(E ll )':x l =2x 5 -4y 12 

(X 2 )' : A = A 

(^i, 2 )': x 5 + 2t l -l = 0 



2t x + 



—2 x 4 

3 



_i_ x 5 

4 



Consider h := 2 and go back to step 4. 



-1 
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4. (F 1J )':-2* 4 -3y 11 =0 
( F 1 J )’: 2»,+^- 1 = 0 

5. Yes. 

6 . - 

7. a = 1 and i 0 = 2 , j 0 = 5 

(A) x 5 = -2t t + 1 statement (V, ) 

(B) Substituting it, we obtain: 

Xl =-6, l+ 2-4y l2 + ^i 

(^ 1 , 2 ) • x 2 ——2t x -\-\ + y l2 
y : x 2 =3tj + l-l + 2y 12 

(C) Remove the equation (f i 2 ). 

(D) Consider v = 2, / = 2 and go back to step 3. 

3. (E u ): x ] =-6t l -4y n +2+^ 

(^2,2 )• X 2 = — 2 t\ >’|2 + 1 
4 - ( F 2.i): -2x 4 -3y 12 = 0 

5. Yes. 

6 . - 

7. n = 2 and i 0 = 1, j 0 = 4 

8 . 2*1 

9. — 3 = — 2 - (1) — 1 

10. (A) x 4 = ->’ 2I + 1 2 statement (// 2 ) 

(B) After substitution, we obtain: 

(£ 2 ,i)’: x, =-6t, ~4>’ l2 + 2+ 2 ’ 2 

(^2,i)': -J 2 i- 2 C=° 

Consider h := 3, and go back to step 4. 

4 - (^ 2 , 1 ) • — ^21 — 2t 2 = 0 

5. Yes 

6 . - 

7. a = 1 and i 0 = 1, j 0 = 21 (two, one). 

(A) v 2] = -2 1 2 statement (V 2 ). 

(B) After substitution, we obtain: 



(C) Remove the equation (F 2 1 ) . 



(D) Consider v = 3, / = 3 and go back to step 3. 
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3. (E 3] ): Xj = -6 t x - 4 y 12 -2t 2 +2 
(^ 3,2 )• -U — —2-t\ + y |2 + 1 



4. 

5. 



No. The general integer solution of the system is: 



Xj = -6fcj - 4 k 2 - 2 k 3 + 2 , 
x 2 = -2k ] +k 2 +1, 

- 1 , 



x 2 — 3 + 2k 2 



x 5 = -2k ] 



3L 



+ 1 , 



from (£ 31 ); 
from ( E 3 2 ) ; 
from (//,)'; 
from 

from (Vj); 



where (k v k 2 
Stop. 



Note 1. Algorithm 3, 4, and 5 can be applied in the calculation of the integer 
solution of a linear equation. 



Note 2. The algorithms, because of their form, are easily adapted to a computer 
program. 



Note 3. It is up to the reader to decide on which algorithm to use. Good luck! 



REFERENCES 

[1] Smarandache, Florentin - Rezolvarea ecuatiilor si a sistemelor de ecuatii 
liniare in numere intregi - diploma paper, University of Craiova, 1979. 

[2] Smarandache, Florentin - Generalisations et generalites - Edition 
Nouvelle, Fes (Maroc), 1984. 

[3] Smarandache, Florentin - Problems avec et sans .. problemes! Somipress, 
Fes (Maroc), 1983. 

[4] Smarandache, Florentin - General solution proprieties in whole numbers 
for linear equations - Bui. Univ. Brasov, Series C, Mathematics, vol. 
XXIV, pp. 37-39, 1982. 

[5] Smarandache, Florentin - Baze de solutii pentru congruence lineare - Bui. 
Univ. Bra§ov, Series C, Mathematics, vol. XXII, pp. 25-31, 1980, re- 
published in Buletinul §tiintific si Tehnic al Institutului Politehnic “Traian 
Vuia”, Timisoara, Series Mathematics-Physics, tome 26 (40) fascicle 2, 
pp. 13-16, 1981, reviewed in Mathematical Rev. (USA): 83e:10006. 

[6] Smarandache, Florentin - O generalizare a teoremei lui Euler referitoare la 
congruente - Bui. Univ. Brasov, series C, mathematics, vol. XXII, pp. 07- 
12, reviewed in Mathematical Reviews (USA): 84j: 10006. 



56 




[7] Creanga, I., Cazacu, C., Mihut, P., Opait, Gh., Corina Reischer - 
Introcucere in teoria numerelor - Editura Didactica si Pedagogica, 
Bucharest, 1965. 

[8] Cucurezeanu, Ion - Probleme de aritmetica si teoria numerelor, Editura 
Tehnica, Buharest, 1976. 

[9] Ghelfond, A. O. - Rezolvarea ecuatiilor in numere intregi - translation 
from Russian, Editura Tehnica, Bucharest, 1954. 

[10] Golstein, E., Youndin, D. - Problemes particuliers de la programmation 
lineaire - Edition Mir, Moscou, Traduit de russe, 1973. 

[11] Ion, D. Ion, Nita, C. - Elemente de aritmetica cu aplicatii in tehnici de 
calcul, Editura Tehnica, Bucharest, 1978. 

[12] Ion, D. Ion, Radu, K. - Algebra - Editura Didactica §i Pedagogica, 
Bucharest 1970. 

[13] Mordell, L. - Two papers on number theory - Veb Deutscher Verlag der 
Wissenschafen, Berlin, 1972. 

[14] Popovici, C. P. - Aritmetica §i teoria numerelor - Editura Didactica §i 
Pedagogica, Bucharest, 1963. 

[15] Popovici, C. P. - Logica §i teoria numerelor - Editura Didactica §i 
Pedagogica, Bucharest, 1970. 

[16] Popovici, C. P. - Teoria numerelor - lecture course, Editura Didactica §i 
Pedagogica, Bucharest, 1973. 

[17] Rusu, E - Aritmetica si teoria numerelor - Editura Didactica si 
Pedagogica, Bucharest, 1963. 

[18] Rusu, E. - Bazele teoriei numerelor - Editura Tehinca, Bucharest, 1953. 

[19] Sierpinski, W. - Ce slim si ce nu slim despre numerele prime - Editura. 
§tiintifica, Bucharest, 1966. 

[20] Sierpinski, W. - 250 problemes de theorie elementaires des nombres - 
Classiques Hachette, Paris, 1972. 



[Partly published in “Bulet. Univ. Bra§ov”, series C, Vol. XXIV, pp. 37-9, 1982, 
under the title: “General integer solution properties for linear equations”.] 



57 




